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Abstract 



In the framework of the functional renormalization group (FRG) we present a simple truncation scheme for the computation of 
real-time mesonic «-point functions, consistent with the derivative expansion of the effective action. Via analytic continuation on 
the level of the flow equations we perform calculations of mesonic spectral functions in the scalar 0{N) model, which we use as an 
D exploratory example. By investigating the renormalization-scale dependence of the 2-point functions we shed liaht on the nature of 
the sigma meson, whose spectral properties are predominantly of dynamical origin. 
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1. Introduction 

(—! 

'"2^' Real-time observables in strongly-interacting systems of- 
I ten pose major challenges for theoretical calculations. Spec- 
tral functions are key to many such observables. They provide 
^ information on quasi-particle spectra and collective excitations 
of the system. The spectral functions of the energy-momentum 
'tensor moreover allow to extract transport coefficients via the 
Kubo formulas HIS HI 01 and thus to study macroscopic prop- 
lerties of strongly-interacting matter 

Functional renormalization group (FRG) methods IH Q I?! 
have been well developed over the last thirty years. Although 
various formulations exist, their underlying physical principles 
are always the same. The equations of motion are described as 
functional differential equations of a scale-dependent generat- 
ing functional and functional derivatives of it. 

Here we employ the approach pioneered by Wetterich 
|with the scale-dependent effective average action as the cen- 
tral object, which generalizes the usual effective action (see 
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1211 for a general introduction). By integrating 



the Wetterich equation which describes the evolution of the ef- 
fective average action with the RG scale k, one obtains the full 
quantum effective action in the limit A; — > 0. The evolution 
equations for the «-point functions can be obtained from that 
of the effective action via functional derivatives. 

A typical problem one encounters in functional approaches 
is the fact that the flow equations for n-point functions involve 
information on up to (« + 2)-point functions thereby leading to 
an infinite tower of coupled evolution equations. Practical ap- 
plications thus require truncations in order to obtain a closed 
system of equations to solve. One frequently used truncation 
scheme is the derivative expansion based on expanding the ef- 
fective action in terms of derivatives. The leading order in this 
expansion is the so-called local potential approximation (LPA). 
Despite their simplicity, derivative expansions in general and 



the LPA in particular, have been applied successfully to a broad 
range of physical systems and critical phenomena. 

In this letter we present an extension of the derivative ex- 
pansion to obtain real-time 2-point correlation functions con- 
sistent with the underlying truncation for the effective poten- 
tial. Starting from the flow equation for 2-point functions we 
employ LPA vertices to obtained a closed system of flow equa- 
tions involving as the only input the scale-dependent effective 
potential. A consistent truncation scheme for the calculation of 
the Euclidean 2-point correlators at finite external momentum 
is, however, only the first step towards a calculation of real-time 
quantities such as spectral functions or transport coefficients. 

The crucial second step, a common challenge in Euclidean 
approaches to thermal quantum field theory, is the analytic con- 
tinuation of the external momentum to Minkowski space-time. 
At present, real-time correlation functions are usually either re- 
constructed from their Euclidean analogues using Pade approx- 
imants or by maximum entropy methods. Alternatively, they 
can be calculated directly in Minkowski space-time by an ana- 
lytical continuation at the level of the flow equations iTsL 141. 
In this respect, our approach is similar to that of Ref. lIlSll in 
which more refined truncations in real time were proposed to 
include effects of the back-coupling of a non-trivial propagator 
on the effective potential and the wavefunction renormalization. 
Actual solutions, however, require an Ansatz for the form of the 
propagator and its singularity structure. Our approach, on the 
other hand, does not rely on assuming a certain spectral shape, 
but it ignores back-coupling effects at present. 

In the following, we concentrate on the calculation of spec- 
tral functions within the 0{N) model. The analysis of the 0{N) 
scalar model, which is often used as a chiral effective model for 
QCD (A^ - 4), has a long history. Spectral functions in 0{N) 
models have been calculated for example in optimized perturba- 
tion theory | 12 , 18 , 3 by taking into account the cr ^ n+n 



process. The critical exponents have been evaluated in 111711 . Al 
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though quantum fluctuations are included via a resummed loop 
expansion, the critical exponents remain at their mean-field val- 
ues. To overcome this limitation a more suitable resummation 
scheme is required. Such a framework is provided by the FRG, 
where the application of the derivative expansion leads to very 
accurate results, e.g., for critical exponents in Of^V) models in 



21j| but to a 



higher orders of the derivative expansion 
surprising degree already also in the LPA 1122 

Spectral functions have been calculated with the FRG in 
non-relativistic models using different truncation schemes in- 
cluding the "BMW" approximation |24ll25n . vertex expansions 
or derivative expansions |26, 27[ l28ll . In all these cases, how- 
ever, the spectral functions were reconstructed from Euclidean 
2-point correlators via analytic continuation using Fade approx- 
imants. Here we employ a simple LPA for the effective action 
which we then use to solve analytically continued flow equa- 
tions for the 2-point functions at real frequency. For simplicity, 
we restrict ourselves to the zero-temperature case in this work. 
In principle it is possible to extend the proposed approximation 
scheme to finite temperature and finite chemical potential or to 
include fermions lfl3[[l4ll . 



2. Formulation 

2.1. Functional RG 

The functional renormalization group provides a powerful 
non-perturbative tool, especially for the study of critical phe- 
nomena. It describes the evolution of a scale-dependent effec- 
tive average action from the microscopic bare action, spec- 
ified at an ultra-violet (UV) scale k - A, to the corresponding 
full quantum effective action for ^ — > in terms of a functional 
differential equation J?!], 



(1) 



where F^"' denotes the «* functional derivative of F/; with re- 
spect to the fields, and Rt is a suitable regulator function. The 
trace is taken over the internal degrees of freedom and momen- 
tum space. Fig. [T] shows a diagrammatic representation of the 



dk^k - 



which is the three-momentum analogue of the LPA-optimized 
regulator ll29ll . The introduction of the regulator suppresses 
the propagation of field modes with momenta smaller than the 
renormalization scale k. Integrating the flow equation ([T) from 
the bare classical action at A: = A down to A: = yields the 
full quantum effective action which includes quantum fluctua- 
tions from all momentum modes. Taking n functional deriva- 
tives of the flow equation ([TJ one obtains the flow equations for 
the «-point correlation functions which, in turn, contain up to 
(« -I- 2)-point correlation functions. 

2.2. Flow equation for the effective action 

The derivative expansion corresponds to an expansion of the 
effective action in terms of gradients 
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(3) 



where the effective potential Uk and the wavefuction renormal- 
ization factors Z, are parametrized in terms of the 0{N) in- 
variant (p^ - (pjcp' and an explicit symmetry breaking is in- 
duced by the ccr term giving rise to a finite pion mass. To 

(2) 

Study the 0{N) symmetry breaking one decomposes F^ for 
constant background 0, into its transverse (F;r) and longitudi- 
nal (Fo-) components in field space 

ifj^(^,p) = F£(0^p)(., - ^) .Fl^(0^,)^. (4) 

Inserting the Ansatz Q into the flow equation ([TJ and evaluat- 
ing it for a constant field configuration yields the correspond- 
ing flow equation for the effective potential. In the simplest 
case of the LPA which corresponds to setting Z\ k{(l>) - 1 and 
^i,k{ip) - for / > 1 one is left with solving the flow equation 
for the effective potential which reads 



^ = i/(i) + (A^_ 1)1/(1) 
dk I -' ^ 'l" 



(5) 



Here we have defined the loop functions 7*"^ as 



/S"^ = Tr„ 



dkRkiq) 



(6) 



Figure 1 : Diagrammatic representation of tlie flow equation for the effective 
action of tlie 0(N) model. Dashed lines represent scale-dependent boson prop- 
agators and gray filled circles coiTespond to regulator insertions duRk- 

flow equation ([T]). Note that although the flow equation has a 
one-loop structure, it is an exact equation, involving full (scale- 
and field-dependent) propagators. 

The precise form of the regulator function Rk is not fixed but 
limited only by certain general requirements [jsl. A convenient 
choice for our purposes is the sharp regulator 

Rk^(k^- t)B{k^ - q\ (2) 



where the trace runs over the loop-momentum q. For the regu- 
lator choice in Eq. (|2]i, /J'' and if'' are given by 



where E„ = VFTItT', = yjk- + 2U' + W'4r and U' = 
■Mr^- The curvatures 2U' and 2U' + 4C/" at the minimum of the 
effective potential yield the screening masses squared of pion 
{nf"~) and sigma meson (nf"~), respectively. They are defined 
as the spacelike limits p — > of the pion and sigma 2-point 
functions, i.e., the inverse of the corresponding susceptibilities. 
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2.3. Flow equation for the 2-point functions 

Applying two functional derivatives to the flow equation 
Eq. ([T]i, we obtain the exact flow equations for 2-point corre- 
lation functions. Their diagrammatic forms for the pion and the 
sigma in the 0{N) model are shown in Fig.|2] These flow equa- 
tions contain scale-dependent 3- and 4-point functions. In order 
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Figure 2: Diagrammatic representation of the flow equation for a mesonic 2- 
point function in the 0(N) model. Dashed lines represent scale-dependent bo- 
son propagators and gray filled circles correspond to regulator insertions diiRk ■ 
Black filled circles and gray filled squares denote the scale- dependent 3- and 4- 
point vertex functions. 

to close the infinite tower of functional equations for the «-point 
functions, truncations are required. A possible scheme has been 
proposed in Refs. 124 12511 (BMW). The idea is that, because of 
the insertion of the cutoff function, the external momentum de- 
pendence of the scale-dependent 3- and 4-point functions has a 
weaker effect than their dependence on the momenta of internal 
lines containing the loop-momentum. One thus expands the 3- 
and 4-point functions in their external momenta in this scheme. 

Here we choose a simpler truncation which is a natural ex- 
tension of the derivative expansion. For this we replace the 3- 
and 4-point functions by their corresponding scale dependent 
but, at the leading order derivative expansion, momentum inde- 
pendent forms as obtained from the flow equation for the effec- 
tive potential, i.e.. 
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We then obtain the flow equations for the pion/sigma 2-point 
functions as follows: 

dtrf^ = {j^, + j„^){4u"4,f 

- UfAU'\N 4- 1) - ^^if'iAU" + W"(/>-) 

J^Al2U"cp + 8C/"'f')2 + J„„i4U"<pf(N - 1) (9) 
i4-'(12C/" +48t/">- + 16C/"'>'*) 
^I^^\N - l)(4-U" + W"(p\ 



with loop functions Jtj defined as 



Jip)ij = 
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where the trace runs over the momentum q. Note that only the 
symmetric components of J are needed in Eq. (|9]l. In LPA, for 
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Figure 3: Diagrammatic representation of the contributions to the flow of the 2- 
point functions. Full gray circles correspond to insertions of dtRii. Black circles 
and gray squares indicate the scale-dependent 3- and 4- point vertex functions, 
here obtained from {S). 



vanishing spatial external momentum p, and with the regulator 
in Eq. these are obtained as 



JipQ)ij + J{pQ)ji = 

(Ei + Ejf(E] + EiEj + E]) + (E^ + E])pl 



(11) 
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AE]E][pl + {Ei + EjY) 



In Fig. [3l we show diagrammatically the contributions to the 
pion/sigma 2-point functions. Although they look like 1-loop 
contributions, the internal lines represent full scale-dependent 
propagators which we evaluate for the constant field value at 
the minimum of the effective potential. These propagators cor- 
respond to resummations of perturbative diagrams. For the pion 
2-point function the process ;r — > ;r + cr is included and for the 
sigma 2-point function the processes cr cr+cr and cr n+n. 
Other contributions, such as cr —> cr -i- tt, would break the resid- 
ual 0{N - 1) symmetry of the broken phase and are excluded. 

In the following we discuss the solutions to Eqs. (|9]l at zero 
temperature. Note that then, no term proportional to 1 /[p^ -i- 
(Ea- - £;r)^] arises because Landau damping, i.e., the absorption 
of an on-shell particle from the heat bath, does not occur 

We reemphasize that the correlation functions at vanishing 
external momentum in our truncation are consistent with the 
effective action since the flow equations satisfy 



dkrfl(p^0)^2dkU',, 



dkTflip = 0) = 25,C/;. + AdkU'^cf . 



(12) 
(13) 



This means that the equivalence between the static screening 
masses from the spacelike p ^ Q limit of the 2-point func- 
tions and from the curvatures of the effective potential at its 
minimum is manifest in our truncation. Moreover, note that 
Eq. (fT2l) automatically ensures the existence of a dynamical 
Nambu-Goldstone boson in the chiral limit. Thus our trunca- 
tion scheme is also "symmetry conserving." 

We solve the flow equation for the real and imaginary parts 
of the retarded 2-point correlators. These are obtained from 
Eq. ^ via the analytic continuation 

rg^w) = f'fipo = -i(co + ie), /? = 0) , for ; = ;r, cr (14) 

which is taken before the evaluation of the flow equations. 
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The spectral function is proportional to the imaginary part 
of the retarded propagator or more explicitly, for a> > 0, 



Pi(a)) = - - 



Imr*->*(w) 



'^(Rerf'«M)%(lmrP«(a;)f 



/ = 7r,cr, (15) 



which is understood to be evaluated in the limit e — > 0^. In the 
numerical evaluations we have to keep a small but finite positive 
imaginary part in the external energy, however 

3. Numerical results 

3.1. Numerical method 

In order to solve the flow equations numerically, we employ 
the Taylor-expansion method. We expand the flow equations (|5]l 
and (|9]i by expanding Uk and the real or imaginary part of F^^. 
(/ = cr, n) around the scale dependent minimum: 

K 

n=0 
L 

L 

Imr<^)«(a;) = ;^<,(a.)(02-^^)«. 



The flow equations for the effective potential and the 2-point 
functions then translate into flow equations for the expansion 
coefficients a„^k, ^J, i^{<^) and cj^ ^(w) in the usual way. Choos- 
ing L = K - \ ensures that the consistency condition (fTZb be- 
tween the pion 2-point function and the effective potential is 
maintained exactly also in a finite Taylor expansion for the nu- 
merical implementation. The flow equation for the sigma cor- 
relator involves up to four derivatives and thus requires K > A. 
We will employ K = 5 and L = 4 in the following. 

The numerical results presented in this article were obtained 
using the Taylor-expansion method as described above, which 
by construction does not take into account the full field-depen- 
dence of the effective potential. Therefore we checked our re- 
sults against the so-called grid method B22fl . which involves dis- 
cretizing the effective potential on a set of grid points in field 
space and solving the corresponding flow equations at fixed val- 
ues in field space, from which the full field-dependent effective 
potential can be reconstructed. The scale-dependent derivatives 
of the effective potential, evaluated at its minimum, can then be 
used as input for the corresponding flow equations (|9} for the 
2-point functions. The disadvantage of the grid method is that 
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3.65 


0.014 
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425.0 


-0.34 


3.40 


0.002 


93.1 


16.4 


299.8 



Table 1 : Parameter sets for a UV cutoff A = 500 MeV coiTesponding to two 
different pion masses. The physical parameters, /;j and the meson masses, are 
given in MeV. 



one incurs considerably larger numerical costs for reaching low 
IR cutoff values than with the Taylor expansions. Fortunately, 
we were able to verify that the results from the grid method 
generally reproduce the Taylor results for k > 0.05 A very well. 

The parameter sets used here for the 0(4) model with UV 
potential of the form Uk=A - a(p^ + b(p'^ are listed in Tab.[T] 

3.2. Pion and sigma meson spectral functions 

In Fig. |4] and Fig. |5] we show the real and imaginary parts 
of the retarded 2-point functions F^ and F* , respectively. These 
were obtained by solving the flow equations separated into 
real and imaginary parts for a small but finite imaginary part 
e ~ 0.1 MeV in the external frequency w + ie as a function its 
real part co. For the pion the imaginary part in Fig. |4] remains 




200 300 
m [MeV] 

Figure 4: Real and imaginary parts of the retarded pion propgator as a function 
of the time-like external momentum w. 

negligible, whereas the real part behaves as demonstrated in 
earlier investigations lfl4[[l3ll . The zero of the real part, which 
corresponds to the pion pole mass, is found at 135.1 MeV as 
compared to a screening mass of 137.2 MeV. This once more re- 
flects the difference between screening and pole masses which 
is a small effect of only a few percent, however, in the purely 
bosonic model lll4lll3[l . At finite density, these values should be 
compared to the pion mass defined via the onset of pion conden- 
sation when coupling the model to an isospin chemical potential 
ifioilTill . The pole mass from our present truncation for the 2- 
point function is considerably closer to the exact mass from the 
Bose-Einstein condensation transition at zero temperature, in 
particular, when fermions are included as first observed in the 
Quark-Meson-Diquark model for two-color QCD 113 1 . 

For the sigma meson, shown in Fig. |5] the imaginary part 
remains negligible but only up to the 2-pion threshold at w « 
2m„ X 275 MeV and tends to finite negative values thereafter. 
The two-pion emission threshold also shows up as a kink in the 
real part. The zero of the real part is found at 350.6 MeV and 
we will argue that this value might be a good estimate for the 
sigma mass even in cases where no clear maximum in the corre- 
sponding spectral function is visible. For comparison, the value 
of the sigma screening mass extracted from the curvature of the 
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Figure 5: Real and imaginary paits of the retai'ded sigma propagator as a func- 
tion of tlie time-like external momentum oj. 



effective potential results to be 425.0 MeV for the parameter set 
with nearly physical pion mass in Tab.[T] 

Fig. |6] displays the pion and sigma spectral functions as 
calculated from Eq. ( fTSl l. The pion spectral function shows a 
sharp peak at 135 MeV to be identified with the pion mass. 
The sigma spectral function exhibits the 2-pion threshold as 
a characteristic increase in the spectral density followed by a 
broad maximum at about 312 MeV above this threshold. Con- 
sidering a case in which the imaginary part Im F*^* were inde- 
pendent of the (time-like) external momentum so that the only 
momentum-dependence was in the real part Re r^^\ it is obvi- 
ous from Eq. (fTST i that the maximum of the spectral function 
would then occur at the zero of the real part with a width as 
determined by the constant imaginary part. This is the case for 
the pion 2-point function where the peak in the spectral function 
and zero of the 2-point function coincide as they must since the 
imaginary part practically vanishes. The same coincidence be- 
tween the sigma mass defined via the maximum of the spectral 
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Figure 6: Spectral functions for the pion (p^) and the sigma meson (p,,-). 



function and the zero of the corresponding real part is of course 
no longer valid for a momentum-dependent imaginary part as 
obtained here, although the zero of the real part may still pro- 
vide an estimate for the sigma mass. For the current parameter 
set the zero of the real part, however, overestimates the sigma 
mass defined via the maximum of the spectral function by 14%. 
This illustrates the difficulties of a mass assignment for a broad 
resonance with a strongly momentum-dependent width. In this 
case the imaginary part does not vanish at the zero of the real 
part. Hence no zero is found in the complex plane of the phys- 
ical sheet, consistent with the Kallen-Lehmann representation, 
which allows only poles on the real axis. One expects a zero on 
the unphysical sheet, however, as illustrated in lll6il3lll . 

Except for a significantly lower sigma mass in our calcu- 
lations, the zero-temperature spectral functions are in qualita- 
tive agreement with perturbation theory results 11161 Il7[ llSl 11911 . 
All calculations show a sharp peak in the pion spectral function 
which determines the pion mass and a sigma resonance peak 
beyond the two pion threshold in the sigma spectral function. 
This agreement is reassuring since perturbation theory around 
the mean-field vacuum was found to work well outside the crit- 
ical region. 
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Figure 7: Comparison of pion (p,r) and sigma meson (p,^) spectral functions for 
parameter sets corresponding to two different pion masses. 

As a consistency check we have also calculated the pion and 
sigma spectral functions closer to the chiral limit using the sec- 
ond parameter set in Tab. [1] corresponding to a pion pole mass 
of 16.0 MeV. The comparison with the physical pion mass is 
shown Fig.|7] When approaching the chiral limit, one expects 
the spectral weight of the pion pole in its correlator to increase 
more and more as this pole moves closer to the oj - Q axis 
where it eventually accumulates the full spectral weight in the 
chiral limit. Consequently, the spectral sum rule then implies 
that all other contributions to the spectral function should de- 
crease with decreasing pion mass. Both these trends are seen 
in Fig. |7] where we extended the frequency range of Fig. |6]to 
include the cr - n threshold in the pion spectral function. To 
explicitly verify the non-renormalization of the pion field in the 
chiral limit requires a more careful analysis, however, proba- 
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bly best done with polar coordinates in field space to correctly 
disentangle the Goldstone bosons from the radial mode 1132 1. 
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3.3. Spectral functions at intermediate scales 

As described in the previous section, we find a resonance 
sigma at A: = 0. It is illustrative to also consider the spectral 
functions at intermediate scales and thus their evolution with 
the RG scale k. In Fig.[8l we show the scale dependence of the 
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Figure 8: Zeros of the real part of tlie sigma 2-point function, the sigma screen- 
mg mass ( V2(7' +402{///) and jhe 2 -pion threshold (2 ^fk^^2U') as a func- 
tion of the RG scale k. 

zeros of the real part of the sigma's 2-point function together 
with its screening mass, defined via as well as 

the 2-pion threshold at 2 ^k- + 2U'. 

During the ^-evolution two distinct zeros are found in the 
real part. We start by considering the first zero (dashed blue 
line) which represents the pole corresponding to a renormal 
ized sigma field. At the ultraviolet scale, the correlator F, 
agrees with the classical one and the zero thus coincides with 
the screening mass, or, in other words, with the bare sigma field. 
This zero gets renormalized by quantum fluctuations and starts 
to deviate from the sigma screening mass at intermediate scales. 

However, it completely disappears at log(A:/A) 2.2 when it 

reaches the 2-pion threshold. 

It turns out that a second zero emerges (solid red line) at 
intermediate scales. This zero is dynamically generated via the 
cr ^ n + n process, and it is therefore not present at the ultra- 
violet scale. The second zero arises above the scale dependent 
cr TT + TT threshold, oj ~ 2E„, decreases with the RG scale and 

eventually becomes constant at log(^/A) 2.2. Thus we find 

that, while the first zero disappears at some scale, the second 
zero survives and corresponds to the zero of the real part shown 
in Fig.|5] 

Fig. E] displays the imaginary part of the sigma correlator, 
evaluated at the first zero (dashed green line), second zero (solid 
red line) and the curvature mass (dotted blue line), respectively. 
At large scales {k ~ A), both the first zero and the screen- 
ing mass are below the scale-dependent 2-pion threshold and 
the corresponding imaginary parts vanish. The finite imaginary 
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Figure 9: Imaginaiy part of the retarded sigma propagator evaluated at the 
scale-dependent 1 st zero (dashed green line), 2nd zero (solid red line) and sigma 
screening mass (dotted blue line). 



part in the numerical calculation is related to the finite imag- 
inary external momentum employed to calculate the spectral 
function, see Eq. ( fT4b . Below log(^/A) ~ -1, the screening 
mass crosses the threshold and the corresponding imaginary 
part of the propagator becomes finite. This behavior should be 
compared to the scale-dependence of the imaginary part in lllSll . 
which corresponds in some sense to an expansion of the 2-point 
function around the screening mass, albeit in a more refined 
truncation scheme. The second zero, which corresponds to the 
zero in the real part that survives until A: = 0, arises above the 
threshold and hence the imaginary part assumes a finite value. 

It decreases with k and becomes constant below log(A:/ A) 2. 

These observations at intermediate scales indicate that the 
final zero in the real part of the sigma 2-point function does 
not emerge from the original single-particle contribution of the 
sigma at the cutoff scale by renormalization effects, but that it 
is predominantly due to the cr ^ jr + tt process. This is in 
contrast to the pion 2-point function whose single-particle con- 
tribution does flow continuously from the ultraviolet to the in- 
frared and thus corresponds to the renormalized original pion 
mass. Such considerations question the reliability of truncation 
schemes based on an expansion around a single scale-dependent 
pole 11511 . It seems that the inclusion of the full external mo- 
mentum dependence of the sigma 2-point function or at least 
an expansion around a more complicated singularity structure 
is required. 

In order to contrast these conclusions for the sigma with 
the corresponding results for the pion correlator, we show in 
Fig. Ho] the analogue of Fig. [8] the zeros of the real part of 
the pion 2-point function. As before, we find two such zeros 
at intermediate scales, the first one (dashed green line) corre- 
sponding to the renormalized single-particle contribution from 
the pion screening mass at the UV scale, and the second (solid 
red line) corresponding to a dynamically generated pion via the 
TT n + a process. In this case the zero at A: = is connected 
to the first zero, however, whereas the second zero vanishes 
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Figure 10: Zeros of the real part of the pion 2-point function and the pion 
screening mass ( V2C/') as a function of the RG scale k. 



around log(A:/A) 1.9. This underlines the fact that the zero 

in the pion 2-point function originates from the ultraviolet pion 
mass by finite renormalization. 



We stress that it is possible to extend our approximation 
scheme to a thermal medium, and to include fermions. We are 
particularly interested in extending the framework to finite tem- 
perature and to investigate the medium-modified spectral func- 
tions. One should be able to observe the asymptotic restoration 
of chiral symmetry also directly in the spectral functions. In the 
region around the phase boundary one might expect to observe 
qualitative differences between the FRG and results based on 
perturbation theory. 
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4. Summary and Conclusions 

In the present work we have extended the truncation scheme 
of Refs. lfl3i[l4ll for the functional renormalization group (FRG) 
equations of 2-point functions. This scheme is consistent with 
the truncation for the effective potential and hence "symmetry 
conserving" by construction. In contrast to earlier approaches 
in the literature, we have employed an analytic continuation on 
the level of the flow equations which allows for a direct com- 
putation of retarded 2-point functions at real frequencies. The 
numerical solutions of the corresponding flow equations yield 
realistic spectral functions of collective excitations, as we have 
demonstrated for the pion and the sigma meson in the (9(4) 
model. Their shape was found to be in good qualitative agree- 
ment with general expectations and previous studies. 

To gain further insight into the origin of the pion and the 
sigma meson, we have evaluated the sigma and pion correla- 
tors at intermediate scales k. It turns out that the pole in the 
pion propagator flows continuously from the initial screening 
mass at the ultraviolet cutoff" to the physical single-particle pole 
for A: — > by the renormalization eff'ects. Since there is no 
renormalization of the pion in the chiral limit this ensures the 
existence of a zero-mode as required by Goldstone's theorem. 
In stark contrast, the sigma pole at the infrared scale originates 
from a dynamical process. We have found two distinct trajecto- 
ries in the singularity structure of the propagator at intermediate 
scales; one corresponding to the renormalized initial mass and 
the other to a dynamically generated complex pole on the un- 
physical sheet via the cr ^ n+n coupling process. In particular, 
the sigma pole at ^ = belongs to the branch which is dynami- 
cally generated. The other branch, connected to the bare sigma 
mass, disappears at some intermediate scale. 
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